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We scrutinize the reliability of covariant and Gram-Schmidt Lyapunov vectors for capturing hydrodynamic 
Lyapunov modes (HLMs) in one-dimensional Hamiltonian lattices. We show that, in contrast with previous 
claims, HLMs do exist for any energy density, so that strong chaos is not essential for the appearance of genuine 
(covariant) HLMs. In contrast, Gram-Schmidt Lyapunov vectors lead to misleading results concerning the 
existence of HLMs in the case of weak chaos. 
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Introduction.- The existence of hydrodynamic Lyapunov 
modes (HLMs) in spatially extended dynamical systems has 
recently attracted a great deal of attention. HLMs are collec- 
tive long-wavelength perturbations associated with the small- 
est positive Lyapunov exponents (LEs). Theymay appear in 
hard-core and soft-core potential systems |ll|-|9[], as well as in 
coupled-map lattices 1 13- 121, of either Hamiltonian or dis- 
sipative type. HLMs are believed to be connected to a sys- 
tem's macroscopic properties and to encode valuable informa- 
tion for understandiing universal features of high-dimensional 
nonlinear systems OTS I3]- Nowadays, we lack a com- 
plete understanding of how or where these modes may show 
up, although, it is generally believed that conservation laws 
and translational invariance are essential for HLMs to ap- 
pear |3, 5, lo jisll . In the case of Hamiltonian lattices it 
has also been argued ifTii ITsIl that the system must be well 
above the so-called strong stochasticity threshold (SST) iflil - 
1811 in order to show significant HLMs. This stochasticity 
threshold separates weak and strong dynamical chaos: Be- 
low the SST (low energy density) the relaxation time to en- 
ergy equipartition grows as a stretched exponential as the en- 
ergy density decreases |i3-18|, whereas above the SST (for 
high energy densities) the relaxation time needed to reach en- 
ergy equipartition is independent of the energy density. Now, 
for the Fermi-Pasta-Ulam (FPU) | l^l and other Hamilto- 
nian lattice models, it has been claimed lfl4l [Tsll that strong 
chaos is essential for the existence of significant HLMs. This 
relation could suggest a connection between the HLMs and 
the ergodic problem, which is closely related to the dynam- 
ical foundations of statistical mechanics in high-dimensional 
nonintegrable Hamiltonian systems. However, this is indeed a 
puzzling situation since, in principle, no matter whether or not 
energy equipartition is reached within the observation time, 
the system is still chaotic with a finite density of non-zero 
LEs. 

In this work we solve this question by demonstrating that 
covariant HLMs associated with nearly zero LEs do actually 
exist in high-dimensional nonintegrable Hamiltonian systems 
for any energy density. We also show that the significance 
measures of these covariant modes are completely equivalent 
both above and below the SST. In contrast, if non covariant 



Gram-Schmidt (GS) vectors are used, as in previous studies, 
the existence of HLMs and their significance do depend on 
the scalar product convention, resulting in artifacts that are 
not intrinsic to the system under study. 

Most studies in the existing literature dealing with the prob- 
lem of computing HLMs relied on the backward (also called 
GS) Lyapunov vectors (BLVs) to obtain the phase-space di- 
rections that correspond to the smallest positive LEs. BLVs 
are obtained by the standard GS orthonormalization procedure 
to obtain the Lyapunov spectrum 12 ll 12211 . BLVs are known 
to have many important issues and artifacts l23l 12411 that ul- 
timately render them useless for many purposes. First and 
foremost, they are forced to form an orthogonal set; this is not 
a minor point because it leads to vectors that are not covariant 
with the dynamics. Therefore, when left to evolve freely with 
the tangent space dynamics, they will not map to themselves 
but will (exponentially) rapidly collapse in the direction of the 
main Lyapunov vector (LV). An even more important issue 
for our purposes is the fact that BLVs depend explicitly on the 
scalar product needed to define mutual orthogonalization. 

We employ here the so-called covariant (or characteris- 
tic 1 251) Lyapunov vectors (CLVs) ll26ll . This set of vectors 
reflects the bona-fide covariant tangent space directions and 
provides a genuine Oseledec splitting of phase space (one di- 
rection is unambiguously associated with each nondegenerate 
LE). A comparison of the significance of HLMs obtained from 
BLVs and CLVs in coupled map lattices was recently made by 
Yang and Radons 12811 . but it did not reveal important differ- 
ences. In contrast, we show here that strong differences appear 
between BLVs and CLVs below the SST. 

Models.- The reference Hamiltonian for the one- 
dimensional lattice models we are considering can be 

written as H = YliLi [Pi/C^^^i) + - Qi)]' where 

N is the system size and V{x) is the nearest-neighbor 
interaction potential. {mi,qi,pi}fL-^ are the dimensionless 
mass, displacement, and momentum of the ith particle, each 
of unit mass rrii = 1; periodic boundary conditions are 
assumed ((7„^ J = q-^). The models hereafter considered are 
(i) the FPU /3-model, with potential V{x) = + x*/A, 

and (ii) the XY model with V{x) = 1 — cos(x). The energy 
density e = E/N, E being the total energy, is the control 
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parameter for the system dynamics. The integration of the 
2N equations of motion as well as the computation of the 
associated CLVs are performed employing a computationally 
efficient numerical algorithm we presented in Ref. ll29ll . 
In particular, CLVs are obtained from the intersection of 
subspaces embedded by GS (backward) and by so-called 
forward Lyapunov vectors, accordin g to the formulas by 
Wolfe and Samelson 12711 (see also 12411 ). Our numerical 
implementation allows us to explore the low-energy region 
wherein exceedingly long transients are required. 

Analysis.- The significance of HLMs is usually an- 
alyzed through the Lyapunov vector fluctuation density, 
which is a dynamical quantity defined as u^°'\q,t) = 
^ where Qi is the position coordinate 

of the ith particle and i^g^-"'' is the coordinate compo- 
nent of the ath LV [(5q(") {t) ; 5p(") (t)] for the ith par- 
ticle. In order to detect the existence of wavelike or 
spatially extended modes, we apply a Fourier transform, 
M(")(fc,<) = J{dq/2Tr)eyip{ikq)u'^°'\q,t), and then com- 
pute the stationary structure factor of u^°'\q^t), S^^^k) = 
limt_i.oo(u^"-'(''^jO^''"''(~^7 0)' where (■ • • ) stands for time 
average. In order to compare the relative weight of the spec- 
trum maximum with respect to the background power we nor- 
mahze by the area below the curve and define S^"'^ (fc) = 
S^") (k) I Efe -S^"^ (fc)- Typically the spectra for LVs associ- 
ated with near-zero positive LEs will exhibit a maximum at 
some fcmax(a)- 

For a quantitative determination of the existence of HLMs, 
two complementary measures can be used ifTll [Til [Tsi [30ll : 
^^"Hfcmax) and the spectral entropy iJ^"). Here ^^^H^max) 
denotes the height of the ath LV stationary structure factor at 
its maximum, whereas the spectral entropy iJ^"' is defined as 
= - Xife 5(")(fc) In S*(")(/c) and measures how the nor- 
malized spectral power S^^^k) spreads among wave num- 
bers k for a given Lyapunov index a. Lower values of the 
entropy indicate that most of the spectral power is concen- 
trated around fewer wavelengths. Therefore, a small value of 
and a corresponding large value of 5^"-' (fcmax) at some 
index a indicate the existence of a sharp, well-defined peak 
of the structure factor for the ath LV, indicating that a partic- 
ular spatial wavelength is favored. HLMs, if present, would 
correspond to extended macroscopic modes with a long wave- 
length, fcniax(a) — > 0, associated with the smallest positive 
LEs, A„ ^ 0+. 

Results.- Our first result concerns the dependence of the 
significance measure of HLMs on the metric used when con- 
structing the set of backward LVs. Since some scalar prod- 
uct always needs to be defined (for Gram-Schmidt orthonor- 
malization) -and this is to a great extent an arbitrary choice- 
meaningful conclusions can only be obtained if the signifi- 
cance measures are independent of the chosen scalar prod- 
uct. In order to emphasize the impact of the scalar prod- 
uct when searching for HLMs, we use the scalar product, 

[5q(") (t);5p(") (t)] M2JVX2JV [5q'") {t);Sp^"'> {t)f, with 
the matrix (M)i.j = Sij + rj Si<Nj<N, rj being an arbi- 
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FIG. 1: (Color online) Above the SST: (a) Normalized ^("'(fcmax) 
and (b) spectral entropy _ff vs. Lyapunov index a corresponding 
to the BLVs (lines) and CLVs (symbols) for different values of the 
perturbation parameter rj of the scalar product (see text). Both plots 
correspond to the FPU model with A'^ — 128 and e = 10. 



trary constant. This corresponds to the perturbed Euclidean 
norm ^/(Jcff^l + ^]) + (6p)^ . We computed the BLVs and 
CLVs for different choices of i] and the results for the mea- 
sures S'"(fcmax) and ij'"-* for the FPU model in the strongly 
chaotic regime for a high energy density e = 10 are presented 
in Fig. [T] For both sets of vectors it is clear that S'"(fcniax) 
presents a maximum and ij'"^ a minimum for some index 
ttmax- The curve in Fig. [Tta) for the BLV with 77 = corre- 
sponds to the mode reported in Ref. |14|] as a truly HLM for 
the FPU system. However, note that amax ^ 102 is still far 
from 126 by 20%, which is the Lyapunov index value coiTe- 
sponding to the smallest positive LE. This deviation was at- 
tributed in Ref. 11411 to fluctuations of finite-time LEs. We 
claim that this mode is not a good HLM. As clearly seen in 
Fig.dta), BLV curves exhibit an a,„ax that progressively shifts 
as 7/ increases, demonstrating that the mode is not intrinsic and 
its position strongly depends on the employed scalar product. 
In contrast, if CLVs are used instead, the positions of the ex- 
trema, as well as the significance measures 5^"-'(fci„ax) and 
are independent of the i] value used, as confirmed by 
Fig. [T] Also, the positions of both extrema are located at LE 
index a,nax ~ 126, exactly where the smallest positive LE 
is located. These results indicate that only CLVs can generi- 
cally detect the existence of intrinsic (scalar-product indepen- 
dent) HLMs. This becomes even clearer if we compute the 
significance measures ^<"'(fcmax) and for energy den- 
sities well below the SST, which is around Sc ~ 0.20 for the 
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FIG. 2: (Color online) Below the SST: (a) S("'(fcn,ax) and (b) spec- 
tral entropy _ff vs. the LV index a for the BLVs (open circles) and 
CLVs (filled circles). Both plots correspond to the FPU model with 
TV — 128, energy density e — 0.01, and 7; = 0. 



FPU Ill4ll (note also an earlier estimate of Ec ~ 0.12 Bill ). 
In this weakly chaotic regime it has been reported that HLMs 
fail to exist HIH. In Fig.|2]we plot 5(")(fcmax) and i?^") 
for an extremely low energy density e — 0.01, where the sys- 
tem behaves almost as a harmonic oscillator chain. As can be 
readily seen, if CLVs are used, the maximum of ^^"•'(/cinax) 
corresponds to the minimum of at the Lyapunov index 
Qmax ~ 120. Thus covariant HLMs exist even well below the 
SST. In contrast, BLVs fail to detect any significant HLMs for 
e = 0.01 (actually, for any energy density below Ec), which 
again indicates that BLVs are unreliable objects whereby to 
detect HLMs. 

We can now quantify the significance of the HLMs as a 
function of the energy density by plotting in Fig.[3]the extreme 
values ^niax and iJmin indicated in Fig.[T] The energy depen- 
dence of these extreme values also shows striking differences 
when BLVs or CLVs are used. Again, different scalar prod- 
ucts lead to different curves if BLVs are used. GS HLMs seem 
to appear only above some energy density. Indeed, a crossover 
toward an increasing significance of HLMs is observed as the 
energy density increases. However, as shown in Fig. [51 the 
curves for BLVs display a crossover that shifts toward lower 
e values as ?/ increases. Thus, the relation of this crossover 



point with the SST Sc, first conjectured in Ref. Ill4ll . turns out 
to be an artifact that arises from the use of BLVs, since its po- 
sition explicitly depends on the employed scalar product. In 
contrast, the maximum power S'niax(e) and minimum spectral 
entropy i/min(£) remain constant for the full range of energy 
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FIG. 3: (Color online) (a) S'max and (b) i^min vs. energy density 
£ for the BLVs using different scalar products (see text): 77 = 
(open circles), 77 — 0.5 (squares), = 1, (diamonds) and rj — 1.5 
(triangles). Data points computed by means of the CLVs are shown 
with filled circles, and vertical dot-dashed lines indicate the position 
of the SST according to Ref. flUl . 



densities if CLVs are used to identify HLMs. Figure [3] sum- 
marizes our main result: covariant HLMs do exist and their 
significance, in qualitative and quantitative terms, is the same 
for the whole range of energy densities studied, irrespective 
of whether we are above or below the SST. 

Next we briefly discuss our results for the XY model con- 
sisting of a one-dimensional lattice of interacting rotors. For 
large enough energy densities, individual particles become 
nearly independent and the system is close to a collection of 
free rotors, whereas for small energy densities the system is in 
a near-harmonic regime, analogous to the corresponding one 
in the FPU model. We have studied the existence of HLMs 
in the XY model and obtained results that are qualitatively 
similar to those in the FPU case. In the high energy density 
regime the interaction potential is not strong enough to main- 
tain the needed coupling between neighboring rotors; thus 
spatially extended perturbations cannot be supported by the 
system dynamics and no HLMs were observed in this regime. 
In the opposite limit of extremely low energy densities we 
found no HLMs if BLVs are used, in agreement with an ear- 
lier study I15II : however, covariant HLMs do exist even for 
a low energy density of e = 0.1, as can be appreciated in 
Fig.ffl The numerical values of the significance measures de- 
pend again on the scalar product if BLVs are used, while no 
effect is seen for CLVs, as expected from our discussion of the 
FPU case. This confirms our finding that strong chaos is not 



4 



0.2 



1 0.15 



0.1 - 



0.05 
3.5 
3.4 
3.3 
'"'3.2 
3.1 
3.0, 



(a) 



O-OBLVs 
• -•CLVs 



€1. 



(b) ■ 

:g- e-0 -G- €)-e-o- ja o-eH3- -o-e-^ 







16 



32 

a 



48 



62 



FIG. 4: (Color online) (a) Maximum power spectrum S'™' (fcmax) 
and (b) spectral entropy //'"^ for the XY model with N — 64, 
£ = 0.1, and 77 = 0, well below the SST. 



essential for the appearance of genuine covariant HLMs. 

Conclusions.- We have studied HLMs in Hamiltonian lat- 
tices and compared covariant with backward LVs. There are 
strong theoretical arguments which suggest that HLMs are as- 
sociated with the existence of conservation laws and transla- 
tional invariance fs, IJ, 13] of the equations of motion. These 
exact symmetries are thus expected to be satisfied irrespective 
of the energy density. We have shown that covariant HLMs 
indeed exist for any energy density so that strong chaos is 
not an essential requirement for the appearance of genuine 
HLMs. Actually, the significance of covariant HLMs remains 
constant both above and below the SST, which seems to in- 
dicate that the mechanisms held responsible of the existence 
of both phenomena should be entirely different. We have also 
demonstrated that BLVs lead to misleading results concerning 
the existence of good HLMs and to significance measures that 
display an undesired dependence on the scalar product. 

The spatial structure of the covariant HLMs herein stud- 
ied indicates that these extended modes have a scale-invariant 
structure, which contrasts with the truly wavelike spatially pe- 
riodic form characteristic of hard-core systems llH9fl. It re- 
mains to be understood whether the extended modes in these 
two types of systems actually belong to the same kind of phe- 
nomenon, and whether these modes can be related to hydro- 
dynamic properties for arbitrary systems with extended chaos. 
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